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1. INTRODUCTION 

A small electrical network operating in island condition is very sensitive towards any disturbance 
due to the mixed energy demand and sources. The stability condition is achieved when the dynamical 
behavior in the generation side is able to sustain for any level of demand and ensured the dynamical 
frequency network within the permissible level. The unbalance power condition may lead to the severe 
frequency deviation problem [1]—[5]. The investigations on the influence of frequency load control on the 
islandic power system network fed by mixed generation sources are reported in [6]-[10]. The unbalance 
power estimation is highly required before performing the network safety actions to avoid the damage of 
local generator and total load black-out. The effect of power system imbalance towards the renewable energy 
penetration and power electronics application was investigated as in [11]—[14]. These studies show that the 
mixed energy and fast growing of electrical technology and strategy may influence the power flow quality 
and system stability. Hence, monitoring the network condition is crucial due to the unpredictable uncertain 
contingencies. One of the power system parameters which reflects the power system dynamic is the total 
power deficiency. 
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Total power deficiency is defined as the magnitude of the imbalance between load demand and 
generated supplied. This magnitude is considered as variable as the state is basically reflect in the network 
frequency dynamic. The basic to identify the power deficit of the power system network with distributed 
generator dependance is referred to the general swing equation. This approach has commonly been used by 
the load shedding scheme to determine the size of load to be shed upon the disturbance through the rate of 
changed of frequency initial slope as reported in [15]-[20]. All of these works are assumed that the system 
structure has the well-known parameters aspecially the equivalent system inertia which may deviated due to 
the variation of network capacity. Thus, the power deficit determination might be not accurate. 

Inertia constant is the parameter which directly influence the frequency dynamics during the 
transient before the steady-state. The investigation of frequency dynamic and system weight response was 
reported in [21]-[25]. These studies proof that the system inertia is the important parameter to be accounted 
when investigating the caused of poor stability condition. Hence, the design of filter using theoretical 
approach is another solution to enable the power deficit estimation to become more reliable in any 
conditions. This paper utilized a generator model operating in island mode which interconnected with loads. 
An optimal H,, filtering problem was chosen as an approach to design an estimator for this network model. 
The design procedure is referred from the work done by Geromel and Oliveira [26], on linear time-invariant 
continuous time system. Then the designed estimator was verified for nominal and robust estimation during 
uncertain inertia constant parameter case. There are two contributions that has been highlighted in this paper 
which are first, the new approach to estimate the total power deficit for isolated power system network using 
Hœ norm. Second is an approach to estimate the input vector by introducing an additional block function 
denoted by h augmented with the isolated power system network model. 


2. METHOD 

The framework of isolated power system network is identical to hydropower system that consists of 
hydro governor, transient droop compensation, hydro turbine and load which can be described by (1)-(3). The 
input for this model is the electrical load demand P, while the state variables are the frequency w, mechanical 
power Pn and the governor power P,,,. The output of this model is the frequency deviation =< Noted that, all 
the variables are in per unit system. Plus, the model does not consider an automatic generation control (AGC) 
so that the dynamical behaviour of the frequency becomes non-conservative. 


dwu) _ 1 

a TH (Pmpu) a Pecpu)) (1) 
dPm 1 
“at = rp Eor — Pm) (2) 
Pow i ce 
Te = > (Pres — 37 Fav) 3) 


The total power deficiency defines the total power deference between load demand power and 
mechanical power of the generator. In this research, the network dynamical frequency response was observed 
and then augmented with the estimator system to estimate the total power deficit. Consider the isolated power 
system network that is represented in linear-time invariant, continuous time system: 


X(t) = Ax(t) + Bug(t) 
y(t) = Cx (t) + Dwg(t) (4) 
z(t) = C,x(t) 


Where x € R” is the state variable, wg E R™ is the disturbance input; y E€ R” is measured output and z € R* 
is the vector to be estimated. The matrices A and B were obtained through (1)-(3). y is the output state 
denoting the frequency dynamic while z is the output state denoting the expected power deficit to be 
estimated. Noted that the power deficit is originally not the state variable. It is related to the electrical 
changes in demand that can be reached from the input state. To make the linear estimator solution become 
feasible, the power deficit related to the input state must be observable. Hence, it must be converted into the 
state variable. Figure 1 shows the block diagram to formulate the estimator design problem. The low pass 
filter denoted as h was introduced to convert the power deficit input state to become as a variable state and 
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augmented with the generator model. The purpose of this approach is to facilitate the design process and to 
ensure the proper estimator transfer function. 


Power 


Wa she : e 
system Estimator, E —> 
network, G = 


Figure 1. The block configuration for design an estimator 


The full state-space realizations with augmented of system ‘h’ can be written as (5): 


X= Anx + Bawa 


y=C,x+D, (5) 
Pa = CgXx 
with 

0 1/2H 0 0 
ees 0 —1/tr 1/tr 0 
n |=1/RT, 0 —1/T, 0 p 

0 0 0 —1/T, 

—1/2H 
paa Ue ih 6 
n n 0 , ( ) 
Bh 


C,=[1 0 0 0] 

Cg=(0 0 0 Ch] 

Noted that T, is the time constant of an additional low pass filter while the other symbol parameters 
are already mentioned. The problem to be dealt with is to design an estimator to estimate Pj of P4 which 


given by Py = F.y, where F belongs to a linear estimator with minimum state space realization in the form: 


Ê =AR+B(C,x + Dy) 
= AR + BrCyx + BD, (7) 


Py = C8 


The matrix Aç E R™*"f, Br E R™*"r, Cp E RS*"f are to be determined. Connecting (5) to (7), it 
corresponds to the overall augmented system that consist of network model x and estimator model %. 


l= a al kl + [e Joa (8) 


The transfer function from noise input wg to the estimation error e can be written as 
x r i 7 x ; i ‘ . 
Tu (s) := C (sI = A) B = GE — hand the matrices A, B, and C of compatible dimensions are given by: 
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The goal for this formulation is to solve the following design problem: 

Hœ estimation problem: to find a guaranteed estimation performance index y such that 
||GE — h||2, < y produces a feasible upper bound y over the estimator state-space realization in (7). This Ha 
norm can be solved accordance to the lyapunov equation solution in schur compliment form as shown in. 


min y 
Subject to 
AP+PA' B PC 
P>0 ; a -I 0 |]<0 (10) 
E m —yl 


Where: P is positive definite matrix 


2.1. Nominal case 

In nominal case, the isolated power system model parameters are assumed precisely known. The 
upper bound of the H,, norm of Ty(s) := GE — h provided in (10) deserves some comments. For precisely 
known parameters, the objective function ||GE —h||2, <y yields the very small norm below the small 
optimal y > 0 (i.e., exactly equal to zero) and does not introduce any convertism. In this case, the overall 
system matrix structure can be converted into the error system through the decomposition approach. From 
(9), the overall augmented system matrix can be decomposed by pre-post multiply with the invertible 


transformation matrix T and T~1. Choosing T = l | Ary = An — BC,, and Cp = C,, the new matrix 
An, B, and C,, is becomes: 


A =t| 4 0 pot _ fan 0 
n= lec, An — BC = Ay — B,C, 


B, =T Br |- Bn 11 
E |s;D,. |e, + B,D, ( ) 
Cr=[-Cg Cg|T-*=[9 Cy] 

Now the overall augmented system matrix has been decomposed and become two separate systems 


which are augmented generator system and error system. By only taking out the error system 
(the second row) and substitute into (10), the Schur compliment matrix can write as (12): 


(An — BeCn)P + P(An — BpCn) Bn t+ByDn PC, 
P>0 ; Py =I 0 
a | -yl 


<0 (12) 


If this matrix is feasible, then the objective function stated in (10) will be satisfied. However, this 
matrix is hard to solve due to the multiplication of the multiple variables. To this end, the matrix needs to 
convert into linear matrix inequality (LMI) form. If this conversion is accomplished, the solution will become 
feasible. 


au ae p-1 
By partitioned P as P := i? 0 4 and pre and post multiply this matrix to (12). Then the matrix 
will become LMI as (13): 
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T YAn +QCn + An'¥+Cn'Q’ YBn—QDn Cg! 
le 1203 a —I 0 |<0 (13) 


The new matrix variables are given by positive definite matrix Y = P~1 and matrix Q. Hence, the 
estimator design problem is equivalent to the programming problem expressed in terms of LMI (13). The 
estimator matrix is defined by Ay = A, — BfCn ; By = -Y-1Q; Cy = Cy. 


2.2. Parameter uncertainty case 

In this case, the isolated power system model has contaminated with the parameter uncertainty. The 
uncertain parameter that has been considered in this case is the inertia constant H. This is the only parameter 
that significantly affect the slope of dynamical frequency and can be clearly seen in (1). Note that the 
parameter uncertainty is confined to a given polytope satisfying M = )f_, 4;M; for some A; > 0 such that 


An B 
A, +++/, =1 and M := ie ea Only A, and B, matrix of the augmented system is assumed to be 
n n 


unknown. 

The guaranteed estimation performance is accordance to the objective function such that 
sup||GE — h||2, < y and produced a feasible estimator which minimize the y over the estimator state-space 
realization in (7). The matrix in (10) can be written as (14): 


AP+PÄ B PË 
P >00; m -I 0|<0;i=1,..N (14) 
E m —yl 


Partitioning the matrix P := f al satisfying the inequality (14) and multiply to the left by 
=" 
y 0 
furthermore multiplying the constrain P > 0 to the left by j’ and to the right by 7, the above inequalities is 
equivalent to: 


j':=diag[j’, 1, 1] and to the right by j and 7: vik Introducing the new variable as explained in [24], 


Z Z 

k v1 > 0; 

ZAn, + An;Z ZAn; + An;Y +C,F'+Q' ZB; C3 —G' 
E YAn; + FC, + An;Y + CF’ YB,, + FD, Cy <0 (15) 
E E —I 0 
E E E —yl 

i=1,..N 


Hence, the estimator design problem is equivalent to the following programming problem expressed 
in terms of LMI on the variable positive definite matrix Z = Z’ and Y = Y’ and matrix Q, G and F. The filter 
matrix is defined by Ar = —Y-*Q(I —Y~*Z)*; Bp = -Y"*F;, C; = GU -Y*Z)7?. 


3. RESULTS AND DISCUSSION 

In nominal case, the inertia constant parameter value is assumed known with H =5 and an additional 
of transfer function h = 1/0.001s + 1. The estimation performance index y such that ||GE — hl|2, < y was 
solved with resulted upper bound y = 8.0067e — 6. The sudden load demand changed at time 3 seconds has 
triggered the turbine system and made the network frequency drop and deviated to the new operating 
condition. Figure 2 shows the estimated total power deficit in time response while Table 1 shows the 
tabulated data related to the estimation error performance at steady-state response through the integral 
absolute error (IAE), integral square error (ISE) and root means square error (RMSE). The result shows that 
the estimator able to estimate the total power deficit accurately with small estimation error. 
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Figure 2. Estimated total power deficit using Hə norm with 0.2 per unit sudden load demand changed 


Table 1. The estimation error performance using H,, norm method under nominal case 
Inertia constant value IAE ISE RMSE 
5 0.0017 _2.75E-06 _2.13E-04 


For the parameter uncertainty case, the inertia constant parameter value H is assumed known but 
with uncertain minimum and maximum values which range 4.5 to 5.5 and an additional of transfer function 
h = 1/0.001s + 1. Under the 0.2 per unit of sudden load demand changed at time 3 seconds, the robust Ho 
estimation was guaranteed accordance to the objective function sup||GE — h||2, < y such that the resulted 
the upper bound y = 0.0459. Figure 3 shows the estimated total power deficit in time response when the 
system contaminated with uncertain inertia constant value. 
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Figure 3. Estimated total power deficit using Ha norm with 0.2 per unit sudden load demand changed and 
uncertain inertia constant values 


The response shows that the H,, estimator is fragile when the inertia constant parameter changed. 
There is offset between the estimated and the reference at the transient and steady-state response. Table 2 
shows the error performance at steady-state response through the IAE, ISE and RMSE. The results proof that 
the error yields the similar performance for each inertia constant value but has about 4.5% of different 
between estimated and the reference total power deficit value. 


Table 2. The estimation error performance at steady state using Ha norm method under robust case 


Inertia constant value IAE ISE RMSE _ Estimated total power deficit 
4.5 0.4412 0.0039 0.0088 
5 0.4425 0.0039 0.0088 0.1911 
5.5 0.4443 0.0039 _ 0.0089 


Bulletin of Electr Eng & Inf, Vol. 12, No. 5, October 2023: 3153-3160 


Bulletin of Electr Eng & Inf ISSN: 2302-9285 O 3159 


4. CONCLUSION 

The total power deficit for isolated power system network was successfully estimated by using the 
Hœ norm method. The model was simulated via MATLAB with Simulink toolbox. The simulation result is 
used to validate the effectiveness of designed estimator and the error performance analysis was focused on 
steady-state response for two cases which are nominal and parameter uncertainty case. The objective function 
was guaranteed with subject to the LMI constrain to obtain the minimum upper bound of the Hœ norm worst 
magnitude in singular value. The results show that the Hœ norm estimator has an ability to estimate the total 
power deficit with very small estimation error in nominal case but not in parameter uncertainty case when the 
network has perturbed by uncertain inertia constant parameter. 
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